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Abstract. 

This article gives necessary and sufficient conditions for the formation of trapped sur- 
faces in spherically symmetric initial data defined on a closed manifold. Such trapped 
surfaces surround a region in which there occurs an enhancement of matter over the aver- 
age. The conditions are posed directly in terms of physical variables and show that what 
one needs is a relatively large amount of excess matter confined to a small volume. The 
expansion of the universe and an outward flow of matter oppose the formation of trapped 
surfaces; an inward flow of matter helps. The model can be regarded as a Friedmann- 
Lemaitre- Walker cosmology with localized spherical inhomogeneities . We show that the 
total excess mass cannot be too large. 
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I Introduction 

We intend to study the geometry of initial data for the Einstein equations coupled 
with a matter field so as to investigate the presence (or absence) of trapped surfaces. A 
trapped surface is a 2-surface which exists at a particular instant of time and has the 
property that all the outgoing light-rays from it converge [1,2]. One expects that outgoing 
light-rays diverge, so one can immediately deduce that the gravitational field in the vicinity 
of a trapped surface must be very strong so as to prevent the light-rays from expanding. 
The presence of a trapped surface indicates that the spacetime is undergoing gravitational 
collapse. One of the singularity theorems of general relativity states that if one has a 
trapped surface (and if the enclosed volume is finite) then there must be a singularity to 
the future (e.g., [2]). If one accepts Cosmic Censorship, that singularities are hidden, then 
the presence of a trapped surface is a sign that a black hole is in the process of forming. 

Initial data for the gravitational field consists of four objects [g a b, K ab , p and J a ], 
where g ab is the three-dimensional positive definite metric of a 3-hypersurface E; E is to 
be regarded as a spacelike slice through the 4-manifold. K ah is a symmetric 3-tensor which 
is the extrinsic curvature of E as an embedded surface, p is the energy density and J a is 
the momentum density of the matter. These data are not independent, they must satisfy 
the constraint equations 

^R[g] - K ab K ab + (K\) 2 = 16np , (1) 
V a K ah - V b K a a = -8nJ b , (2) 

where ( 3 'i?[<7] is the scalar curvature of E and V is the 3-covariant derivative compatible 
with g ab . 

2 



The trace of the extrinsic curvature (trK = g ab K ab ) is equal to the (positive) time rate 
of change of the three-volume, -^(d 3 V) = trK dV. It is very common to place conditions 
on trK so as to select a preferential slicing ("extrinsic time"). In asymptotically flat 
space-times a standard choice is the maximal slicing condition, trK = 0. In cosmological 
models(e.g. [3]), the analogous condition is to choose trK = constant. 

A trapped surface is defined as a compact two-dimensional (smooth) spacelike surface 
S having the property that the expansion 9 of outgoing future-directed null geodesies which 
are orthogonal to S is everywhere negative on S. If S is regarded as a sub manifold of E, 
then 9 can be expressed in terms of g ab and K ab by 

9 = V a n a - K ab n a n b + g ab K ab , (3) 

where n a is the outward unit normal to S in E. 

We have recently derived simple necessary and sufficient conditions for the appearance 
of trapped surfaces in asymptotically flat initial data sets [4] and in open universes [5] , with 
the additional assumption that the initial data are spherically symmetric. In this article 
we will derive equivalent necessary and sufficient conditions for trapped surfaces in a closed 
universe. We continue to demand spherical symmetry and we impose the condition that the 
data satisfy trK = constant. The problem we wish to consider looks like a spherical lump 
superimposed on a standard Friedmann-Lemaitre- Walker (FLW hereafter) background, in 
the fact that we assume that the density becomes constant outside some finite subset 
and the matter current vanishes outside the same subset. We address also the situation 
where a closed universe is filled with many spherically symmetric bumps that have a small 
domain of influence so that in between them geometry coincides with the homogeneous 
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and isotropic FLW geometry . 

The standard time-slice through a closed FLW cosmology is defined by a three- 
manifold which is just the standard round three-sphere, i.e., the line-element can be written 

as 

ds 2 = a 2 [dr 2 + sin 2 rdn 2 } , (4) 

where a is a constant; K ah is pure trace, i.e., K ab = Hg ab , where if is a constant, the 
Hubble constant; the energy density is a constant, po, and the matter current density is 
zero. The scalar curvature of the manifold is given by 

«>R = £ , (5) 
so therefore the Hamiltonian constraint reduces to 

— + 6H 2 = 16npo . (6) 

or 

The momentum constraint is trivial. 

The normal n a to the two-surfaces of constant radius is given by n a = (a -1 , 0, 0) and 
its divergence is given by 

_ „ / . 2 N _i^sin 2 r 2 

V a n = (asm r) ; = — cotr. (7) 

dr a 

Therefore the expansion 9 of such surfaces is given by 

e= 2 C otr + 2aH (g) 

If the universe is expanding (H > 0), then trapped surfaces can only occur on the 'other' 
side of the equator, i.e., r > where cotr is negative. Since cotr becomes unboundedly 
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large and negative, we will always find trapped surfaces. These trapped surfaces indicate 
the existence of the final 'big crunch'. This property should hold for any spherically sym- 
metric universe and so finding trapped surfaces near the 'south pole' is not very interesting. 
Thus we will be interested in finding trapped surfaces in the 'northern' hemisphere, far 
away from the 'south pole'. These, we expect, should indicate the onset of some local 
gravitational collapse, superimposed on the overall expansion and subsequent collapse of 
the whole universe. 

The initial data [g a b, K ab , p, J a ] we wish to consider are defined on a three- manifold E, 
which has topology S s . We assume trK = g a bK ab = 3H = constant holds everywhere on 
the 3-manifold. We further assume that the metric, the extrinsic curvature, p, and J a are 
spherically symmetric. Because of the spherical symmetry we know that the three-metric 
is conformally flat and so can write the line-element in isotropic coordinates as 



where a is a positive constant and is a positive conformal factor that depends only on r. 
We assume that p and J a are arbitrary (other than satisfying some positivity conditions 
one may wish to impose). We do not impose any equation of state on the matter. We 
need not do so as we restrict attention to a single instant of time, we make no use of the 
Einstein evolution equations. We further assume that p = po (a constant) and J a = 
outside some subregion of compact support. We choose the constant a in eqn.(9) to satisfy 



ds 2 = a 2 ^[dr 2 + sin 2 rdil 2 } , 



(9) 



-z + &H 2 = 167TP0 



(10) 



a 



so as emphasise the link with the FLW cosmology. 
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II The extrinsic curvature in spherically symmetric closed universes 

We assume that the extrinsic curvature is spherically symmetric. This means that the 
general form of K ab is 

K ah = K(r)n a n b + b(r)g ab , (11) 

where n a is the unit vector in the radial direction, and K(r) and b(r) are two scalar 
functions. Since we know trK = 3H, we can rewrite (11) as 

ab 

R ab = Rg ab + R ^ ^a n b _ _ (12) 

o 

The momentum constraint, Eq. (2), reads 

2 lK ' + (6^ + 3 ^)K] = - S -4, (13) 



3a 2 (p 4 <p sinr a<p 2 

where we assume J a = Jn a : with n a = (a _1 ^ _2 ,0,0) as the unit vector in the radial 
direction, J = J(r) is a scalar and the primes represent derivatives with respect to r. Eq. 
(13) can be slightly simplified to 

K' + (6^ + 3 — )K = -UTva^J . (14) 
(p smr 

Now let us make our standard assumption, that J is zero outside some finite subset. 
The momentum constraint (Eq. (14)) reduces to 

, . d>' cosr N 
K'+(6^- + 3- — )K = 0, 15 
sin r 

in the exterior region. It is easy to solve this equation to give 

K = C((j) 6 sin 3 r)~\ (16) 
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where C is some constant. This is a well-known result, in disguise. K generates a 
divergence-free, tracefree (TT) tensor. The spherically symmetric closed manifold is con- 
formally flat, and TT tensors are conformally covariant. There is a unique spherically 
symmetric TT tensor on flat space, and the tensor K(n a n b — g ab /3), with K given by Eq. 
(16), is exactly this tensor, transformed in the correct way. 

The constant C is essentially the integral of the matter current. Eq. (14) can be 
rewritten as 



One important consequence of (18) is that if the matter is at rest, i.e., J = 0, then 
we can immediately deduce that K = and the extrinsic curvature is pure trace. Another 
consequence is that is that if we have a globally spherically symmetric model with the 
matter current confined to a region near the north pole, the integral in Eq. (18) when 
we let R approach tt must go to zero. Otherwise, the extrinsic curvature term in the 
Hamiltonian constraint will have a term that diverges like sin -6 r as r approaches re. In 
other words, if J vanishes outside a finite region (and we have global spherical symmetry) 
so also will K. 

Ill The Hamiltonian constraint in spherically symmetric closed universes 

We now can put Eqs. (9) and (12) together and evaluate the Hamiltonian constraint. 
We can write the scalar curvature as = 6a~ 2 <j)~ 4 — 8</>~ 5 V 2 </>, where the laplacian is 



[Kcjf sin 3 r]' = -127ra0 8 sin 3 rJ . 



(17) 



This can be integrated from r = to r = R to give 




(18) 
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with respect to the background metric (Eq. (4)). We also have that (trK) 2 — K ab K ab = 
QH 2 — \K 2 . Thus the hamiltonian constraint, Eq. (1), reads 
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(19) 



a 
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Rearranging, and using Eq. (10), we get 



V 2 </> = ^(<f> - </> 5 ) - [±-K 2 + 2n(p - po)]^ . 



(20) 



This is our key equation and much of our effort will be devoted to analysing it. 

The other important quantity we need to consider is the expansion 9 (Eq. (3)) of the 
surfaces of constant r. The normal to such surfaces is given by n a = (a -1 </>~ 2 , 0, 0), and 
V ' a n a is given by 



On using Eq. (12) we get -K ab n a n b + g ab K ab = 2H - \K. Hence we get (from (3)) 



The aim of this article is to take the expression (22) for the expansivity in terms of the 
'potentials' (f> and K, and use equations (18) and (20) to replace these quantities with p 

— * 

and J, more physical objects. Two other spherically symmetric cases have already been 
successfully dealt with in this fashion, the case where the three-geometry is asymptotically 
flat and the case where the system asymptotically approaches a flat Friedmann cosmology 



V a n a = (a < /) 6 sin 2 r)- 1 ( ( /) 4 sin 2 r)' = 2(a^ 3 sinr) -1 ^' sinr + 0cosr) . 



(21) 



9 = 2(o(/) 3 sinr) 1 (2(p' sinr + (p cos r) + 2H - -K . 



(22) 



[4,5]. 
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IV. A sufficient condition for trapped surfaces when J = 0. 

We assume that we are given spherically symmetric initial data [g a b, K ab , p, J] (which 
may be confined to a finite region). We assume that the trace of the extrinsic curvature 
is constant. We further specialize by assuming that the matter-current is instantaneously 
zero. As we have shown in Section II, this allows us to conclude that the extrinsic curvature 
is pure trace, so that the function K(r) = 0. We also assume that p assumes a constant 
value po outside a region of compact support. The Hamiltonian constraint now reduces to 
a simplified version of Eq. (20) 

V 2 ^= i ^(^-0 5 )-27rA^ 5 , (23) 

where Ap = p — p . 

We now wish to get some relationship between our sources (in this case Ap) and 
the expansion 9 of any given spherical surface, with coordinate radius r = Rq in the 
northern hemisphere. The trick is to take Eq. (23), multiply it by 4> and integrate it in 
the background metric, over the volume V enclosed by the surface at Rq. This gives 

/ (/}V 2 (/)dV = ^- [ (c/} 2 - cj) 6 )dV - 2tt / Ap<f> 6 dV . (24) 
Jv 4a Jv Jv 

The left-hand-side of (24) can now be written as 

/ 4>V 2 4>dV = I (fN a (t)-dS a - [ (V(j)) 2 dV 
Jv Jr Jv 

i-Ro 

= Anacpcp'sm 2 r\ Ro - Ana / ((f)') 2 sin 2 rdr . (25) 

Jo 

Now consider the expression (22) for 9, dropping the term in K. This can be rearranged 
to give 

(7ra 2 </> 4 sin 2 r) (9 — 2H) = Anacfxf)' sin 2 r + 2na4> 2 sin r cos r . (26) 
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We should recognise that 4na 2 (f> 4 sin 2 r = A is the proper area of a sphere of coordinate 
radius r in the physical space. We now merge (26) with (25) to give 

f A r R ° 

/ (j)V 2 (j)dV = -(0\ Ro - 2H) - Ana / ((f)') 2 sin 2 rdr - 2na(j) 2 smR cosR . (27) 
Jv 4 J 

We can simplify the right-hand-side of (24) when we recognise that 

f (/) 6 dV = V, (28) 
Jv 

where V is the proper volume enclosed by the surface, and 

Ap^dV = AM , (29) 



v 



where AM is the mass excess as measured in the physical space. We now combine (24) 
and (27) to give 



A f R ° AH 2>V 

-e\ Ro =na / [3(f) 2 sin 2 r+4(^) 2 sin 2 r+2(</> 2 sin r cos r)']dr^ — — -2ttAM . (30) 

4 Jo ' 2 4a 

Consider the integrand of the integral on the right- hand-side of (30). It is 

I = 4> 2 sin 2 r + 4(</>') 2 sin 2 r + 400' sin r cos r + 2(f) 2 cos 2 r (31) 
= 2(f) 2 - (f> 2 sin 2 r + A(f)' sin r(4>' sin r + (f> cos r) . (32) 

Let us now make two (as yet unjustified) assumptions. We will assume 

(f)' < and 0'sinr + (pcosr > (33) 

on the support of the matter. In later sections we will show that these conditions follow 
naturally from the Hamiltonian constraint (20), if we assume that the excess mass Ap is 
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positive. Almost identical conditions have been derived and used in [4] and [5]. We now 
can deduce 

I < 2<P 2 . (34) 

Thus we obtain the inequality we want 

A r R ° A H 31/ 

-9\ Ri) < 2ixa J cj) 2 dr + — - — - 2nAM . (35) 

We can identify a 4> 2 dr = L as the proper radius of the sphere of coordinate radius 
Ro. Thus we can write inequality (35) as 

We have proven the following: 

Lemma 1: Assuming conditions (33), if we can find in our physical data a spherical 
surface satisfying 

r AH W 
AM > L + — — — 37 

4tt 8yra 2 

then that surface is trapped. 

If we were interested in minimal surfaces rather than in trapped surfaces we can prove: 

Lemma 2: Under the conditions stated in this section, if we can find in our physical 
data a spherical surface satisfying 

1 : A»)fl>L-^ 5 , (38) 



167T 87ra 2 

where A^R is the excess integrated scalar curvature, then the manifold must contain a 
minimal surface. 
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V. A necessary condition for trapped surfaces when J = 

It is possible to find a necessary condition for trapped surfaces under much weaker 
conditions. We need make no assumptions such as (33). Let us return to equation (30) 
and consider again the integrand (31) 

2 2 / / \ 2 * 2 / 22 

I = <j) sin r + 4(</> ) sin r + 400 sin r cos r + 2<fi cos r 
= <p 2 + [(pcosr + 20' sin r] 2 (39) 
> 2 . (40) 

It is interesting to note that we do not impose any conditions on 4> or (/)' to derive this 
inequality. When (40) is substituted back into (30) we get 

A r R ° A H "\V 

-e\ Ro >naJ^ <j> 2 dr+ — -^-2nAM. (41) 

Again we use /(f a<p 2 dr = L, the proper radius of the sphere to simplify (41) to 

A L AH 3V 

Thus we have 

Theorem I: Any spherical surface in the physical data satisfying 

L AH 3V 

AM < — I ^ 43 

2 4tt 8yra 2 V ; 

cannot be trapped. 

Let us stress again that this inequality is valid for any spherical surface in the whole 
universe. As a trivial application, it is compatible with the fact that all the surfaces near 

12 



the south pole are trapped, even with no extra mass. As we approach the south pole, 
in the standard Friedmann model, we get L ~ <27r, V ~ 2n 2 a 3 and A ~ 0. Thus the 
right-hand-side of (43) approaches — ira/A. 

VI. The effect of a matter current on trapped surface formation 

In Sections IV and V we have dealt with the situation where the matter current was 
at rest. In this section we will consider the effects of nonzero J. We will continue to 
assume that both the excess matter and the current density are confined to the northern 
hemisphere and that the trace of the extrinsic curvature is a constant. The hamiltonian 
constraint has the general form as given by (20) 



and we see that the effect of the nonzero current is just to add another positive term to Ap. 
We need to derive an equation similar to (30), but taking into account the extra terms that 
arise in both the definition of the expansion and in the Hamiltonian constraint. Before 
we do so, it is useful to multiply the momentum constraint by the unit radial vector and 
integrate it over a spherical volume in the physical space. We assume that the extrinsic 
curvature has the form given by (12) 



V 2 </> = JL(0 - </> 5 ) - [±-K 2 + 2nAp]<p 5 , 



(44) 
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ab 



K ab = Hg ab + K(r)[n a n b 



3 



(45) 



where H is a constant. We thus wish to consider 




(46) 
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The term in H contributes nothing to the integral so we can ignore it. Thus we can replace 
K ah in the integral with k ab = K(r)[n a n b — g ab /3]. The second term on the left in (46) is 
zero as well. Now we integrate by parts the remaining term on the left to give 

An a n b k ab \ Ro - [ k ab V a n b dV = -8tt / n a J a dV . (47) 
Jv Jv 

This can be further simplified to read 

\kA\ Ro + \ ( KV a n a dV = -8tt / n a J a dV . (48) 

Substituting from Eq. (21) finally gives 

2 d>na C R ° f 

-KA\ Ro + — Kacj) 3 sin r(2<f>' sin r + cos r)dr = -8tt / n a J a dV . (49) 

3 3 Jy 

Let us now return to the Hamiltonian constraint. Repeating the manipulation of Section 

V leads to an equation analogous to (30) 

A f R ° 1 AH AK W 

-e\ Ro =naJ^ [I - -{Ka^sinrf]dr + - |^ - 2ttAM , (50) 

where I is exactly the integrand given in (31). Let us eliminate the AK term from (50) 
by subtracting 1/4 of (49) from it. To simplify the notation somewhat, we will replace 
Kacj) 3 sin r by k and divide by 2tt to give 

A a f R ° 12 AH W 

8^ k = 2~/ lI-zk 2 + -k(2ct>'smr + ( t>cosr)}dr+ — -— ^ - (AM - A J) , (51) 

where we define A J = f y n a J a dV . 

We can write the total integrand in (51) as 

1 

I' = (p 2 (sm 2 r + 7/3 cos 2 r) + 80' sinr(2<^' sinr + </>cosr — k/S) — —(k — 60' sinr — (pcosr) 2 . 

o 

(52) 
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The first term in the integrand I' is easy to handle. It can be written as 

h = </> 2 (7/3-4/3sin 2 r) < ^ 2 , (53) 

and the integral of this term is < 7L/6. The middle term can be written, on using (22) as 

h = 8a^y sin 2 r(0/2 - H) (54) 
= AaOcp'cj) 3 sin 2 r - 8Hci(j) 3 (j)' sin 2 r (55) 
= 4a0<£'0 3 sin 2 r - 2Ha((p 4 sin 2 r)' + AHacf) 4 sin r cos r . (56) 

If we have no trapped surfaces in the interior, we have that 9 > and therefore the first 
term in (56) is obviously negative, on using our standard inequality (33), i.e., (/)' < 0. The 
second term integrates to give — AH/A-k, which cancels the equivalent term in (51). The 
only term that we need to consider carefully is the third term in (56). The integral of this 
term can be estimated as follows. Let us define 

f R 

T = 2Ha 2 / (p 4 sin r cos rdr . (57) 



Jo 

We can show that V < HL 2 . To obtain this we need to use r _1 sin r cos r < 1. This allows 
us to write 

r R 

r < 2H / {arct) 2 ){a(j) 2 dr) . (58) 
Jo 

If we assume </>' < we have that 4> monotonically decreases so we know that, at any point, 
a<t> 2 r < L, and, of course, acj) 2 dr = dL. This immediately gives us the desired result. 

Let us now return to Eq. (52) and assume that there is no trapped surface inside the 
radius Ro we are considering. We can now write 

A 7 3V 

-<, W <-L + HLi-— 2 -(AM-AJ). (59) 
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This now will give us our desired result: 

Lemma 3: Assuming 0' < 0, if we obtain a surface for which 

7 

AM -AJ>-L + HL 2 - (60) 

6 8iva 2 

either the surface is itself trapped or the interior contains trapped surfaces. 

Another estimate can be derived if we are willing to use both conditions in (33). We 
can write the total integrand I' in (51) in a different way. We get 

16 1 

/' = 4> 2 (sin 2 r + 7/3 cos 2 r) H 4>' sinr(</>' sinr + </>cos r) (k — 20' sin r — (f> cos r ) 2 • (61) 

o o 

If we accept (33) the middle term is negative, the last term is obviously negative, and, as 
in (53), the first term is less than |^> 2 . This means that we can replace (59) with 

A 7 AH "\V 

This gives us 

Lemma 4: Assuming conditions (33), if we obtain a surface for which 

7 AM "W 

AM-AJ>-L + --— 2 (63) 

the surface is trapped. 

Note: Lemma 4 is usually better than Lemma 3 because if we assume (j)' < we 
immediately get A < 4itL 2 . 
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VII. The total excess mass must be bounded 

There are many inequalities that one can derive using these techniques, but one, in 
particular, is enlightening i n that it shows that there are very strict bounds on the total 
amount of excess mass that can be placed in finite region, independent of whether this 
matter is inside a horizon or not. Let us return to eqns.(23) and (24) but assume that J 
is nonzero. These can be written as 

/f R ° 3V 
K 2 dV + 2n AM = -Ana^cp' sin 2 r\ Ro +na J (3(p 2 sin 2 r + 4((p') 2 sin 2 r)dr- . (64) 

This can be rewritten as 

J K 2 dV + 2tvAM = -4%a(j) sin r((j)' sin r + 4>cosr)\ Ro + 

f R ° 3V 

2tt acj) 2 sin r cosr \r + ixa I Idr — , (65) 

Jo 4a 



where I is the same integrand as in Eq. (31). The first term is negative (on assuming 
(33)), and the integral is bounded above (from (34)) by 2rcL. Thus we have 

/3V 
K 2 dV + 2nAM < 2-n acj) 2 sinr cos r\ Ro + 2nL- — . (66) 

a<f) 2 sin r = R is the natural areal coordinate, the equivalent of the Schwarzschild coordi- 
nate. Therefore we get 

3V 1 f 

AM < R cosi? + L — -—r - — / K 2 dV. (67) 

8na z 24% J 

Hence the mass of the inhomogeneity inside any given sphere cannot exceed the sum of 
the proper radius and the areal radius. This is a result that supports Einstein's view that 
"matter cannot be concentrated arbitrarily" [6]. 
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VIII. Nonhomogeneous cosmologies. 



In Sections IV, V and VI we have derived necessary and sufficient conditions for the 
appearance of trapped surfaces in spherical cosmologies. The only assumptions we make are 
contained in equation (33). This section, and the following ones will be devoted to showing 
that these conditions can be derived as a consequence of the Hamiltonian constraint (20), 
assuming that the mass excess be positive and localized. 

We can conceive of two very different situations in which spherical symmetry could 
be assumed in a cosmological context. The local situation would be where the spherical 
symmetry holds only on a (small) patch of the data, a spherical galaxy in a universe with 
many such objects. We will discuss this case first. The other case is when the spherical 
symmetry holds globally, we have only one spherical lump in the whole universe. We 
postpone discussion of this situation until Section IX. 

In the last few sections we were concerned with universes containing an isolated spher- 
ically symmetric lump. There may be a large number of such bumps, some of them with 
excess energy but many with deficit energy. Direct integration of the Lichnerowicz equa- 
tion (Eq. (20)) shows that the total energy of all bumps vanishes. If the universe on 
average is homogeneous and isotropic we expect that there exist regions with positive and 
negative excess energy. 

We assume that those regions are separated so that in between them the metric 
approaches the Friedmann-Lemaitre- Walker metric, i.e., the conformal factor <p tends to 
1 outside the perturbed region and the extrinsic curvature vanishes (except, of course, 
for the Hubble constant term). We place the north pole at the center of the given lump. 

18 



Then one might be more specific and say that <p tends to 1 before reaching the coordinate 
distance r = tt/4 or even r = 57r/12. That is not a severe restriction since, as pointed 
above, there might exist many bumps and a zone of influence of any of these may cover 
only a small fraction of the whole universe. 

Keeping this in mind, we prove the following: 

Lemma 5. Assume that we are given an r , < r < 5tt/12 such thatAp is nonneg- 
ative for r e (0, ro) and that 0(ro) > 1 then 

0>1, 0'<O Vre(0,r ) (68) 

Proof. Using spherical symmetry, we may write equation (20) as follows: 

V 2 </> = - 1 2 (sin 2 r0')' = -27rAp0 5 - K 2 /12<p 5 + 3/(4a 2 )</>(l - </> 4 ) . (69) 
a 2 sin r 

One easily checks that 

oj = sup [3/ (4) 0(1 - </> 4 ) : < (j> < 1] = 0.6/5 1/4 = 0.4012442 . 

We will need below an equation related to (66) , 

V 2 x = 1 2 (sin 2 r X ')' = «o • (70) 
a z sin r 

Assume the contrary to the claim: let there exists a solution of (24) such that (j)' > on 
an interval (i?*, i?i). Let us remark that this means that < 1 in a part of this interval 
because if 4> > 1 on the whole interval, then the right-hand-side of (66) is nonpositive. The 
maximum principle guarantees that <p cannot have an interior minimum, which contradicts 
the positivity of </>'. 
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Assume that (p(R*) = %(-R*); by definition d/dr(p(R*) = 0. The equation (67) is 
solved by 

X = -(« /2)r cot r + C, . (71) 

One can find that x is increasing everywhere, so that at the point R* we have d/drx > 
d/drcf). Under the above stated initial conditions, and recognising that V 2 (%-</>) > 0, one 
may conclude that everywhere on the chosen interval 



d(f)/dr < dx/dr, 4> < x ■ 

The change 5(f) = (fi(Ri) — </>(-R*) °f 4> on [Rt-> Ri] 1S estimated from above by the change of 
X, given by x(R) ~ x{R*i) (note that both functions are equal at R* !). But the change of 
X on all intervals {Ji(R* , Ri), n < oo on which <fi is increasing is smaller than the change 
of x on the entire interval (0, 57r/12), which is not greater than a/3. Since at that point 
<p > 1, it means that (/)>! — ao/3. 

Now, let us observe that the last term on the right hand side of (24) is decreasing for 
4> > 1 — «o/3. Define 

ax = sup[30(l - </> 4 )/4 : 1 > > 1 - <V3] , 

according to the above remark, ol\ < «o- Repeating the above described procedure, we 
will obtain a better estimation from below for the function 0, 

Repeating the above procedure infinitely many times, we finally arrive at the desired 
estimation 

(f) > 1 . 
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From Eq. (69) we see that 4> is monotonically decreasing, since the right hand side of 
(69) is nonpositive and we may apply the maximum principle. That ends the proof of the 
lemma. 

Remark. That iteration procedure that has been described above works until one 
meets a fixed point of the map: 

<f> i+1 = 1-^(1 -0*)/4, 

i. e., a value 0* such that 

0* = l-0(l-<^)/ 4 . 

The only solution of the last equation is 0* = 1. 

One situation that is of interest is where one has an isolated lump. This is the 
situation where the deviation from the Friedmann background is localized in the sense 
that the conformal factor equals 1 somewhere in a region outside the lump. Such an 
isolated lump may be formed by having a region of enhanced density surrounded by a 
region with diminished density In such a situation we can prove 

Lemma 6: Given three radii ro,ri,r2 with ro < r\ < ro < 5tt/12 and r\ < it/2, 
such that Ap > for r G (0,r ), Ap < for r G (ro,ri) and Ap = for r G (ri,r 2 ). If 
(p = 1 for a particular r G (ri, r 2 ) and if </> < (sinr) -1 / 2 then > 1, 0' < for r G (0, ri). 

We will postpone the proof of this Lemma until the end of Section X. 

Lemma 6 and Lemma 3 can now be combined to prove 

Theorem II: Given an isolated spherical lump with a central region of positive excess 
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mass and if in that region we find a surface satisfying 



7 "W 
AM — AJ > -L + HL 2 



6 8yra 2 
we must have a trapped surface. 

We will also prove the following. 

Lemma 7.1f Ap is nonnegative and if ((//sinr + <pcosr)\ ro > then 

4>' sin r + (pcosr > 0,Vr G (0, r ) . (72) 

Proof. Notice that 

1 1 

— ~ (0' sinr + </>cosrY = — ^— (sin 2 re/V — d>/a (73) 

a 2 sin r a 2 sin r 

= - 27 rAW> 5 - ( A + 3 ^) - K 2 5 /12 . (74) 

The last equation follows from the constraint equation (24). The quantity <f>' smr + 4>cosr is 
decreasing; but at r$ it is > 0, hence (/)' sinr + cos r is positive everywhere for r G (0, ro). 

IX Spherically symmetric constant scalar curvature manifolds 

Until now we have been dealing with localized spherically symmetric lumps which 
we regarded as being embedded in a Friedmann background. The class of spherically 
symmetric cosmologies contains many other solutions. In particular, we wish to investigate 
solutions which are globally spherically symmetric. The general situation we wish to 
consider consists of a spherically symmetric initia 1 data set which has a lump of excess 
matter surrounded by a constant density background which fills the rest of the universe. 
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Outside the lump, the Hamiltonian constraint guarantees that the scalar curvature is 
constant. The combination of spherical symmetry and constancy of the scalar curvature 
guarantees that that section of the manifold can be regarded as part of a standard round 
S 3 [7] . This is not equivalent to guaranteeing that the conformal factor must be identically 
one in the exterior region. We need to understand the conformal structure of spherically 
symmetric constant scalar curvature manifolds. 

Let us return to the Hamiltonian constraint, Eq. (19) 

V 2 </> = ^ (0 - </> 5 ) - + 2tt(p - p )]^ , (75) 

and let us consider the simplest possible case. Let us assume that the matter is at rest, 
which gives us K = 0, and that the energy density is a constant, p = p . The Hamiltonian 
constraint now simplifies enormously, to give 

V 2 0=^(0-0 5 ). (76) 

Obviously, one solution of (76) is </> = 1, which is the standard solution. This is not the 
only one. There is a complete family of regular solutions of Eq. (76) given by 

i 

(X 2 

0(a) = — r^— — , 2n i , (77) 

(a 2 cos 2 § + sm 2 §)2 

where a is an arbitrary positive constant. 

In Eq. (76) we are trying to find a conformal factor 4> which transforms a manifold 
with constant scalar curvature = 6/a 2 to another manifold with the same constant 
curvature = 6/a 2 . The Yamabe Theorem [8] states that every Riemannian manifold 
can be conformally transformed into on e with constant scalar curvature. The existence 
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part of the Yamabe theorem has bee n shown to be true. Whether the resulting constant 
scalar curvature manifold is unique is an open problem. However, it is known that in the 
case of flat space we do not have uniqueness; there is at least a one parameter family of 
inequivalent manifolds. The conformal factors mapping between them are just the 0( a )'s 
given by Eq. (77). The manifold with metric 

dS 2 = (j)\ a) a 2 [dr 2 + sin 2 rdQ 2 ] (78) 

can be shown to be a round metric by using the coordinate transformation [7] 

L = n — 2 tan -1 (a cot -) . (79) 

It is important that we understand the general structure of the 0( a )'s. It is clear from 
(77) that 0( a ) with a = 1 satisfies <j> = 1 and so represents the identity transformation. 
Each 0( a ), with a ^ 1, equals \ j \fa at r = 0; equals [2a/(a 2 + l)] 1 / 2 , which is always less 
than 1, at r = n/2; and equals -Jo. at r = n. If a < 1 then 0( a ) monotonically decreases, 
if a > 1, 0( a ) monotonically increases. The value of r, call it R\, at which each 0( Q ) passes 
through 1 is given by 

R 

a = tan 2 — - . (80) 

In the range < r < n/2 the 0( a )'s with a > 1 give a smooth set of non-intersecting 
curves with one curve passing through each point in the (r, </>) plane with 0( a ) in the range 
< <j> < 1. At each fixed value of r as a increases, 4> decreases. Each individual curve 
is increasing as r increases but they are all still less than 1 on reaching r = n/2. Beyond 
r = it/2 each curve in turn crosses the <p = 1 line at increasing values of R\ with increasing 
a (as given by (80)) without intersecting any other curve. However, each 0( a ) is climbing 

24 



more rapidly than the curves with lower values of a by the time it crosss the line = 1. 
Above this line the curves proceed to cross because a curve with higher values of a must 
rise above all the curves with lower a by the time it reaches r = it. 

Superimposed on this pattern is a mirror-image set of curves for values of a < 1. 
These start with values greater than 1 at r = 0, the smaller a the greater the value. They 
then proceed to decrease as r increases, crossing one another as the curves with smaller 
values of a fall below the curves with larger values of a. All the crossing is accomplished 
by the time they cross the 0=1 line, and this all happens before one reaches r = tt/2. 
Beyond r = tv/2 all the curves smoothly decrease without crossing until one reaches r = 7r. 
The part of (0, r) space defined by 4> 2 sin r > 1 contains no curves. Through every other 
point two curves, with different values of a, pass. 

Now let us consider one of these constant scalar curvature manifolds, with line-element 
given by (75). When we substitute the explicit form of 0( a ) into Eq. (21) we get 

V„n a = 2(aasinr)- 1 [a 2 cos 2 (r/2) - sin 2 (r/2)] . (81) 

This quantity has only one zero at r = Re given by 

a = tan-^. (82) 

An alternative form is 

1 - a 2 

cos R E = - ~ • (83) 

1 + or 

This value (Re), from (79), corresponds to L = ir/2, i.e., the physical equator of the 
conformally transformed space. V a n° is positive for all r < Re and is increasingly neg- 
ative for r > Re- Therefore, independent of the value of a, all the manifolds share the 
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property of the standard Friedmann slice that they have no trapped surfaces in the 'north- 
ern' hemisphere, but that as one approaches the 'south' pole all the two-spheres become 
trapped. 

The equator, where a = tan^p, also plays another role which will be important in 
our future discussions. Consider one of these curves 0( a ), with a given a, and let us assume 
a > 1. Start on this curve at r = 0, with = 1/y/a < 1. As one moves along this given 
curve, other 0( a )'s, with different values of a, cross it. These curves will have a < 1. These 
begin with curves whose value of a is essentially zero. As r increases, as we move along 
our chosen curve, the value of a of the curves crossing it monotonically increases, but is 
still less than 1 as we pass r = tt/2. Finally, at Ri, given by Eq. (80), with Ri > 7r/2, 
the curv e crosses the line 0=1 (which is, of course, the 0( a ) with a = 1). The chosen 
curve now starts to overtake curves with values of a > 1 but less than our chosen a. This 
continues until one reaches the the equator, r = Re, given by eqns.(79) and (80). Note 
Re > Ri for a > 1. This point also satisfies (f) 2 smR E = 1 so is the boundary of the 
excluded region. Only one curve, our chosen one, passes through this point. From this 
point on , as we move along our curve in a direction of increasing r, our chosen curve is 
now being overtaken by curves 0( a ) with larger and larger a's. Just as the vertical line 
r = is essentially the curve 4> a with a = 0, the vertical line r = it is the </> a with a = oo. 

X Spherically symmetric models with non-constant scalar curvature 

We wish to construct a manifold which is spherically symmetric with a given (non- 
constant) scalar curvature R. We can assume that the metric is of the form (9), and then 
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we can write the equation satisfied by in a form very similar to equation (75), i.e., 

V 2 0= i ^(0-^ 5 )-^A^ 5 , (84) 

where AR = R- 6/a 2 . 

As we mentioned earlier, we wish to consider manifolds with constant density outside 
some compact region. Thus we will wish to consider situations where AR is zero outside 
some finite spherical volume. In the exterior region the equation for reduces to 

V 2 </>=^(</>-</> 5 ). (85) 

We know a complete family of solutions to this equation, the 0( Q )'s of Eq. (77). However, 
it is not immediately obvious that the solution of Eq. (85), which is only valid on part of 
the manifold, must be one of these 0( a )'s. In one case we need regularity both at r = 
and at r = tt. In the other case we need regularity only at, say, r = tt. We know that 
any solution to (84) must have some finite positive value at r = 7r, and that it's first 
derivative must vanish there. We have a 0( a ) with the same value at r = tt and with the 
same first derivative. This would be usually enough to prove that the two functions must 
coincide. We cannot immediately deduce this here because the point r = tt is a singular 
point of the equation because sinr = there. Happily, Eq. (85) i s really only an ordinary 
differential equation because of the spherical symmetry and a treatment exists of such 
singular points in Rendall and Schmidt [9] . Using their Theorem 1 allow us to deduce that 
the two solutions must coincide and that the solution of Eq. (85) must be one of the 0( a )'s 
in the exterior of the support of AR. 

One obvious situation we would like to consider is when there is a step-function in 
the scalar curvature, when the scalar curvature has one (constant) value in part of the 
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manifold and has a different (constant) value on the rest. It is easy to see that the 
conformal factor to achieve this must made up of a pair of the functions 4>{a) defined by 
Eq. (77). More precisely, we would want that the conformal factor equal Ci0( Ql ) in one 
region and C20(a 2 ) m the rest, where C\ and Ci are prescribed constants depending on the 
values of the scalar curvature and ol\ and a<i are adjustable parameters. One now tries to 
match these functions and their first derivatives at some radius and discovers that this can 
not be done, irrespective of the values of Ci, C<i or the coordinate value of the matching 
point.. 

This means that we cannot expect to be able to solve Eq. (84) with a randomly 
chosen AR. Equally, if we get a solution it may well not be unique. Nevertheless, we 
have been able to extract a number of interesting and useful properties of solutions to Eq. 
(84) (assuming one exists!). The situation we are interested in is when AR > and has 
compact support. In the exterior region we know that the solution must be one of the 0( a ) 
's. We will show that this a must satisfy a < 1 if the support of AR is not too large. 

The assumption we will make is that the support of AR lies entirely in the 'northern' 
hemisphere. We do not define this in terms of the background geometry, but rather in 
terms of the physical metric. The exterior solution is defined by a </>( a ) with a = a - We 
assume that the 'equator', the radius Re satisfying 

a = tan R E /2 , (86) 

lies in the exterior zone. 

The analysis is based entirely on the maximum principle as applied to Eq. (84) . The 
solution (p to (84) cuts through a complete family of 0( a )'s. At any point along 0, we know 
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that the two 0( a )'s that pass through that point satisfy 

V 2 ^) = ^(^)-i))- (87) 
Therefore we must have at that point since we assume AR > 

V 2 (0 (a) -0) >0. (88) 

If AR > at this point then by continuity V 2 (</>( a ) — 0) > close to this point. This 
implies, due to the maximum principle, that if 4> approaches a 0( a ) from above it must 
pass through it; it cannot either 'bounce' off it or merge with it. 

We obtain our result by contradiction, so we begin by assuming that «o > 1. Let us 
start at r = 7T and move along the curve <p, the solution to (84). In this region it coincides 
with (f> ao . We continue along this curve through Re (as defined by (85)). At some radius 
r < Re we enter the support of AR. At this point the curve 4> must bend downwards, 
away from (f) ao . But in this region the space underneath (f) ao is filled with all the curves 
0( a ) with a > ao- Our solution curve must continue descending faster than these 0( a )'s. 
As we point out above, it can neither merge with, nor bounce off, any of them. This can 
only lead to disaster. Either <p goes to zero, or hits the r = point with a non-zero angle 
(if (f)' = at r = 0, it would be tangent to one of the 0( a )'s). Thus the original assumption, 
that cto > 1 cannot be true. A similar disaster befalls a curve which starts with ao = 1. 
Therefore we must conclude that a < 1. 

Each point of (r, 4>) space (excepting the region defined by </> 2 sinr > 1) has two curves 
0( a ) passing through it. Thus any curve <f>(r) can be described by giving the pair of a values 
corresponding to each point it passes through. Consider the particular curve we describe 
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in the preceeding paragraph. As it moves away from r = n one of the a's remains fixed 
at «o, the other starts at oo and monotonically decreases, at Re both a's coincide, and 
at values of r < Re the varying a is less than a . Now we enter the support of AR. 
The curve 4>{r) now drops below the curve 4>{a )- One a continues to decrease, the other 
increases above «o- The key question is what happens to this increasing a ? It cannot 
reach a maximum value either at or before reaching r = 0, because this is equivalent to 
being a tangent (from above) to one of the 0( Q )'s, which is forbidden by the maximum 
principle. Thus something bad must happen to the curve in question, it either goes to zero 
at a nonzero value of r, or 4>' > at r = 0. 

We can now deduce somewhat more assuming, as always, that AR > and that it has 
support only on one side of the equator. We now assume cko < 1. Let us follow the solution 
to (84) in from r = n. We start off at a value of </> = (ao) 1 / 2 < 1 and it monotonically 
increases as we move along ao . We will not deviate from this curve until we, at least, 
reach the equator (Re)- By this point 4> ao has risen above 1. The solution cf> starts curling 
downwards relative to (f) ao at some r < Re- However, it cannot curl down too much, 
because if it passes, going downwards, through = 1 it is stuck in the same trap as before. 
We therefore can assume that <fi > 1 in the interior. But , as long as is > 1, we have 
V 2 4> < and this means that <fr cannot have a minimum in the support of AR. Hence we 
must have 4>' < in the interior, with equality only at r = 0. Of course, since «o < we 
have 4>' < in the exterior. 

The result that the conformal factor outside the support of the excess matter must 
coincide with one of the special functions 0( a ) is interesting insofar as it shows that the 
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exterior region is homogeneous and isotropic, uninfluenced by the lump. This is exactly 
analogous with the result derived in Section II that the current potential K vanishes outside 
the support of the current. 

Let us now return to the promised proof of Lemma 6 from Section VIII. Here we 
consider a situation where part of the manifold consists of a region (0, ro) near the north 
pole with AR > 0, a second region (r ,ri) with AR < and that we have standard 
FLW data outside these regions. This means that the conformal factor, the solution to 
(20), equals 1 outside r±. If we follow the conformal factor in from the right it must 
increase as it enters the region of negative Ap. If it went down there would be an interval 
(ri — 5, r\ + 5') in which <fi < 1. But we have from (83), since AR < 0, 

V 2 <p > . 

However <p achieves an interior maximum at r = r\ , which is forbidden. Therefore we must 
have 4> > 1 as we approach r±. 

Having shown that cf> > l,(f>' < near r\, we now wish to show that > 1 in the 
whole interval (ro,ri). Let us assume that (j)' = somewhere in the interval (ro,ri). If 
(j) 2 sin r < 1 in the interval, this means that 4> must turn over in the region that is filled 
with </> a curves. In particular, this means that the 4> curve must tangent from below one 
of these curves, say ao , at some point, say r = r%. At this point we have 

0(r 3 ) = </>a (r 3 ) > 1, 0'(r 3 ) = <f>a (r 3 ) < . (89) 
In a neighbourhood of r 3 we also have 

<P < 0a o • (90) 
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Subtracting (87) from (84) gives 

V 2 (0 ao - 0) = ^[0 Qo - 0* o - + 5 ] + ^AR<j) 5 . (91) 

The right-hand-side of (88) is non-positive, so the minimum principle says that ao — 
cannot achieve a minimum at r^, which contradicts our assumptions. Hence 0' < on the 
whole interval (ro,ri), with equality only in the FLW region. Therefore 0(ro) > 1. Now 
we can use Lemma 5 to finally show > 1, 0' < in the interval (0, ro). 

XI Trapped surfaces in spherical universes. 

In the previous section we have proven that a spherical bump of compact support 
confined to the northern hemisphere produces a monotonically decreasing conformal factor 

0, 

0'<O, 

if its excess energy is non-negative. This corresponds to Lemma 5 of Sec. VIII. We can 
also prove 

0' sin r + cos r > 

Outside the bump must coincide with Q , where a < 1. As is shown in Lemma 7, the 
quantity 0' sinr + 0cosr monotonically decreases. Outside the matter it is equal to 

a -1 3 [(i + cos 2 r)(a 2 - 1) + 2cosr(l + a 2 )] . 

This expression vanishes at a value Ri such that cosi?i = (1 — a)/(l + a), (Eq. (79)), i.e., 
when Q crosses through 1. Therefore inside a region with no n- negative matter density 
the inequality (73) holds. 
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Thus we can prove: 

Theorem III: Given a globally spherically symmetric initial data set which satisfies 
trK = constant with a localized positive lump which is instantaneously at rest. If we can 
find a spherical surface satisfying 

A „, r AH 3V 

AM > L + — 

47T 87ra 2 

then that surface is trapped. 
We can also prove: 

Theorem IV: Given a globally spherically symmetric initial data set which satisfies 
trK = constant with a localized positive lump may be moving. If we can find a spherical 
surface satisfying 

A „ . T 7 T AH 3V 
AM - AJ > -L + 



6 4tt 8yra 2 
then that surface is trapped. 

In addition, one can prove another version of a sufficient condition (Lemma 3) for 

the formation of trapped surfaces by moving inhomogeneities. Indeed, let us return to the 

equation (68). In section VII we estimated V from above by HL 2 . The other estimate 

is that T < 8Ha 2 . This is obtained when one realises that <j> merges with some <f> a , with 

a < 1 in the exterior, and that it must lie below this 4> a in the interior. In turn, we know 

4> a < a -1 / 2 . Thus, if we replace </> 4 in (82) with a~ 2 we will get an upper bound 

rR 

T<AHa 2 a 2 / sinr cos rdr = 2Ha 2 a 2 sin 2 R . (92). 
Jo 

This, as it stands is not a particularly interesting inequality until one realises that we 
are assuming that the excess matter is confined within the equator and that we are not 
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interested in finding trapped surfaces beyond the equator. Therefore, we can replace R 
in (92) by Re as given by (83). It is a straightforward manipulation of trigonometric 
functions to show 

a" 2 sin 2 R E = 4 < 4 , (93) 

(1 + a z y 

which gives the desired inequality. Using (89) and proceeding as in Section VI, we obtain 
the following version of Lemma 3: 

Theorem V: Under the preceeding conditions, if we obtain a surface for which 

either the surface itself is trapped or the interior contains trapped surfaces. 

XI Concluding remarks. 

Let us summarize the whole discussion. In this paper we study spherically symmetric 
closed universes. There are two different situations which are of interest. 

First, we may have a homogeneous and spherically symmetric background geometry 
with many spherical bumps placed in it. Globally, a resulting universe is neither spherically 
symmetric nor homogeneous, but locally, close to a particular lump, there is a spherical 
symmetry with respect to the centre of the lump. Also, very far from the bump, the geome- 
try of the chosen Cauchy slice coincides with the homogeneous and spherically background 
geometry. That case was studied in Section VIII. In Section VIII is formulated a sufficient 
condition (Theorem II) for the formation of trapped surfaces by moving perturbations, 
when the initial momentum of the gravitational field is changed. Obviously, Theorem II 
can be specialized to the case where the excess matter is at rest. Theorem I in Section V 
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gives us a nice necessary condition for the formation of a trapped surface when the matter 
is at rest. Interestingly, we have failed to find an equivalent necessary condition in the case 
where the matter is moving. In Section VII we prove that in a sphere of a fixed radius 
only a finite amount of perturbed energy can be packed. 

Secondly, one might be interested in investigating geometry of a universe that is not 
homogeneous but is globally spherically symmetric. Section IX contains a description 
of spherically symmetric constant scalar curvature models. These solutions are used in 
Section X to derive a number of properties of initial data for nonhomogeneous distributions 
of matter. Using them, we prove in Sec. XI criteria for the formation of trapped surfaces 
in globally spherically symmetric geometries that are generated by a single bump confined 
to only one hemisphere. 

We would like to stress that these results obtained here are correct in the full nonlinear 
theory, nowhere do we draw on a perturbation expansion or a linear approximation. We 
expect that these results are true also for a class of nonspherical perturbations, for instance 
in the class of conformally flat perturbations of the homogeneous closed universes we would 
expect to obtain results analogous to those known in open flat universes [10]. 
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